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Abstract 

A combined eflective model reproducing the equation of state of hadronic 
matter as obtained in recent lattice QCD simulations is presented. 

PACS numbers: 12.38.Gc, 12.38.Mh, 12.40.Ee, 24.85. +p 



1 Introduction 

We are going to construc10 a combined effective model reproducing the equa- 
tion of state of hadronic matter as obtained in recent lattice QCD simulations 
[21 [3] . The model should reproduce basic physical characteristics of processes 
encountered in the dense hadronic matter, from the hot QCD phase through 
the critical temperature region till the lower temperature hadron resonance gas 
phase. It has been shown in [i] that the equation of state derived from that 
time QCD lattice calculation [5] can be reproduced by a simple hadron gas 
resonance model below critical temperature Tc- For higher temperatures the 
model is modified by introducing finite widths of heavy hadrons ^6] with a 
heuristic ansatz for the spectral function which reflects medium modifications 
of hadrons. This Mott-Hagedorn type model is constructed to fit nicely the 
lattice data, also above Tc where it does so because it leaves light hadrons be- 
low a mass threshold of mo = 1 GeV unaffected. The description of the lattice 
data at high temperatures is accidental because the effective number of those 
degrees of freedom approximately coincides with that of quarks and gluons. 
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In order to remove this unphysical aspect of the otherwise appealing model 
one has to extend the spectral broadening also to the light hadrons and thus 
describe their disappearance due to the Mott effect while simultaneously the 
quark and gluon degrees of freedom appear at high temperatures due to chiral 
symmetry restoration and deconfinement. In the present contribution we will 
report first results obtained by introducing a unified treatment of all hadronic 
resonances with a state-dependent width Ti{T) in accordance with the inverse 
hadronic collision time scale from a recent model for chemical freeze-out in a 
resonance gas [?]■ The appearance of quark and gluon degrees of freedom is 
introduced by the Polyakov-loop improved Nambu-Jona-Lasinio (PNJL) model 
[Slin]- The model is further refined by adding perturbative corrections to 0{as) 
for the high-momentum region above the three-momentum cutoff inherent in 
the PNJL model. One obtains eventually a good agreement with lattice QCD 
data, comparable with all important physical characteristics taken into account. 



2 Extended Mott-Hagedorn resonance gas 

We introduce the width F of a resonance in the statistical model through the 
spectral function 

r • TTi 

A{M, m) = N„-, —5 J , (1) 

where Nm is the standard normalization factor. 

The model ansatz for the resonance width F is given by |5] 

where Cr = 10-^ N^n = 2.5, Nt = 6.5 and Th = 165 MeV. 

For simplicity, we assume ns — for the strangeness number density and 
Ub = for the baryon number density. Then the respective chemical potentials 
jjB = and /is = always, so that the temperature is the only significant 
statistical parameter here. 

The energy density of this model can then be cast in the form 

e(r)= 9r£riT]m,)+ ^ I ^d{M^)A{M,m,)e^{T;M), (3) 

mi<mo mi>mo "'o 

where mg = 1 GeV and the energy density per degree of freedom with a mass 
M is 

f <fik Vk^ + A'P 

with the degeneracy gi. For mesons. Si — —1 and for baryons Si — 1. For our 
hadronic gas kept at fixed volume we have the relation 

dP 

e = T-s-P = T- — -P, (5) 
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where P = P{T) and s = s{T) are the pressure and entropy density. 




T [MeV] T [MeV] 

Figure 1: Left: Thermodynamic quantities for the old Mott-Hagedorn Reso- 
nance Gas model [B] . Different line styles correspond to different values for the 
parameter Nm in the range from Nm = 2.5 (dashed line) to Nm = 3.0 (solid 
line). Lattice QCD data are from Ref. [2]. Right: Thermodynamic quantities 
for the new Mott-Hagedorn Resonance Gas where quark-gluon plasma contri- 
butions are described within the PNJL model including as corrections (dashed 
lines). Hadronic resonances are described within the resonance gas with finite 
width, as an implementation of the Mott effect (dash-dotted line). The sum of 
both contributions (solid lines) is shown for the energy density (thick lines) and 
pressure (thin lines) in comparison with the lattice data from [2] . 

In the left panel of Fig.[T]the results for the pressure and energy density of the 
model at this stage are shown. Although providing us with an excellent fit of the 
lattice data, the high-temperature phase of this model is unphysical. Imposing 
that all mesons lighter than mo = 1 GeV are stable provides us with a SB limit at 
high temperatures which mimicks that due to quarks and gluons in the case for 
three flavors [10]. In reality, due to the chiral phase transition at Tc, the quarks 
loose their mass and therefore the threshold of the continuum of quark-antiquark 
scattering states is lowered. The light meson masses, however, remain almost 
unaffected by the increase in the temperature of the system. Consequently, they 
merge the continuum and become unbound - their spectral function changes 
from a delta-function (on-shell bound states) to a Breit-Wigner-type (off-shell, 
resonant scattering states). This phenomenon is the hadronic analogue |11) 
of the Mott- Anderson transition for electrons in solid state physics (insulator- 
conductor transition). It has been first introduced for the hadronic-to-quark- 
matter transition in jl2j . Later, within the NJL model, a microscopic approach 
to the thermodynamics of the Mott dissociation of mesons in quark matter has 
been given in the form of a generalized Beth-Uhlenbeck equation of state [13] , 
see also [H] . 

As a microscopic treatment of the Mott effect for all resonances is presently 
out of reach, we introduce an ansatz for a state-dependent hadron resonance 
width Ti{T) given by the inverse collision time scale recently suggested within 
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an approach to the chemical freeze-out and chiral condensate in a resonance gas 

m 

UT) = r-i,(T) = ^ i^')T{r^)T n,{T) , (6) 
j 

which is based on a binary collision approximation and relaxation time ansatz 
using for the in-medium hadron-hadron cross sections the geometrical Povh- 
Hiifner law [13]. In Eq. (jS]) the coefficient A is a free parameter, nj{T) is the 
partial density of the hadron j and the mean squared radii of hadrons (r|)T 
obtain in the medium a temperature dependence which is governed by the (par- 
tial) restoration of chiral symmetry. For the pion this was quantitatively studied 
within the NJL model and it was shown that close to the Mott transition 
the pion radius is well approximated by 

•iiT) = = ^mn-' . (7) 

Here the Gell-Mann-Oakes-Renner relation has been used and the pion mass 
shall be assumed chirally protected and thus temperature independent. 

For the nucleon, we shall assume the radius to consist of two components, a 
medium independent hard core radius ro and a pion cloud contribution (T) = 
Tq + r'^{T) , where from the vacuum values r-jr — 0.59 fm and rjv = 0.74 fm 
one gets ro = 0.45 fm. A key point of our approach is that the temperature 
dependent hadronic radii shall diverge when hadron dissociation (Mott effect) 
sets in, driven basically by the restoration of chiral symmetry. As a consequence, 
in the vicinity of the chiral restoration temperature all meson radii shall behave 
like that of the pion and all baryon radii like that of the nucleon. 

The resulting energy density and pressure behavior is shown in the right 
panel of Fig. [T] This part of the model we call Mott-Hagedorn-Resonance-Gas 
(MHRG). When all hadrons are gone at T 250 MeV, we are clearly missing 
degrees of freedom! 



3 Quarks, gluons and hadron resonances 

We improve the PNJL model over its standard versions [HI IS] by adding per- 
turbative corrections in 0{as) for the high-momentum region above the three- 
momentum cutoff A. In the second step, the MHRG part is replaced by its 
final form, using the state-dependent spectral function for the description of the 
Mott dissociation of all hadron resonances above the chiral transition. The total 
pressure obtains the form 

P{T) = PmHRg(T) + PpNJL,MF(r) + P2{T) . (8) 

where PuHRGiT) stands for the pressure of the MHRG model, accounting for 
the dissociation of hadrons in hot dense, matter. 

The 0{as) corrections can be split in quark and gluon contributions 

F2 (T) - A""""" (r) + P|'"°" (T) , (9) 
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where p^^^^ stands for the quark contribution and Pf contains the ghost 
and gluon contributions. The total perturbative QCD correction to 0{as) is 

TT TT 

where = J^^rp ■ The corresponding contribution to the energy density 
is given in standard way by Eq. ([5]). 

We will now include an effective description of the dissociation of hadrons 
due to the Mott effect into the hadron resonance gas model by including the 
state dependent hadron resonance width ^ into the definition of the HRG 
pressure 

Pmhrg(T) J -^dM A{M)Tln{l + S^e-^P'+^'Z^) . (11) 

From the pressure as a thermodynamic potential all relevant thermodynamical 
functions can be obtained. Combining the corrected meanfield PNJL model 
for the quark-gluon subsystem with the MHRG description of the hadronic 
resonances we obtain the results shown in the right panel of Fig. [T] where the 
resulting partial contributions in comparison with lattice QCD data from Ref. [2] 
are shown. 

We see that the lattice QCD thermodynamics is in full accordance with a 
hadron resonance gas up to a temperature of ~ 170 MeV which corresponds 
to the pseudocritical temperature of the chiral phase transition. The lattice 
data saturate below the Stefan-Boltzmann limit of an ideal quark-gluon gas at 
high temperatures. The PNJL model, however, attains this limit by construc- 
tion. The deviation is to good accuracy described by perturbative corrections 
to 0{as) which vanish at low temperatures due to an infrared cutoff procedure. 
The transition region 170 < r[MeV] < 250 is described by the MHRG model, 
resulting in a decreasing HRG pressure which vanishes at T ^ 250 MeV. 

We have presented two stages of an effective model description of QCD ther- 
modynamics at finite temperatures which properly accounts for the fact that in 
the QCD transition region it is dominated by a tower of hadronic resonances. To 
this end we have further developed a generalization of the Hagedorn resonance 
gas thermodynamics which includes the finite lifetime of hadronic resonances 
in a hot and dense medium by a model ansatz for a temperature- and mass 
dependent spectral function. 
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